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It is known that the near-field spectrum of the local density of states of electromagnetic field 
above a SiC/air interface displays an intense narrow peak due to the presence of a surface polariton. 
It has been recently shown that this surface wave can be strongly coupled with the surface plasmon 
of graphene in graphene-SiC heterosystems. Here, we explore the interplay between these two 
phenomena and demonstrate that the spectrum of the LDOS in these systems presents two peaks 
whose position depends dramatically both on the distance to the interface and on the chemical 
potential of graphene. This paves the way towards the active control of the LDOS. 



PACS numbers: 73.20.Mf, 07.79. Fc, 44.40.+a, 78.20.-e 

Graphene has attracted during the last years a strong 
interest because of its unique electronic [IH3] and opti- 
cal properties U that could lead to breakthrough tech- 
nologies in the domains of nanoelectronics and nano- 
optics. One of the main features is the linear form of 
graphene electronic dispersion curves which is the sig- 
nature of charge carriers of zero effective mass, travel- 
ing over long distances without losses. More recently, 
the study of graphene optical properties has provoked a 
very particular attention [SHH]- An important property of 
graphene sheets is the presence of strongly confined sur- 
face plasmons on its surface which can be easily tuned by 
modifying the charge density within the lattice of carbon 
atoms [lOT - fl^ . It has been suggested that these tun- 
able surface modes could be used to manage near-field 
radiative heat transfer [131 [IB] or to convert near-field 
energy with graphene-based thermophotovoltaic devices 
|17l[TS] . When studying these surface excitations, the in- 
trinsic quantity containing all the relevant information is 
the local density of states (LDOS). Direct observation of 
surface plasmons on graphene has been reported recently 
[191 [20] using near- field detection techniques. The con- 
nection between LDOS and near-field measurements is 
not straightforward and results from a convolution with 
a complex aperture function. A discussion for different 
setups can be found in refs. pUl - E^ . LDOS due to sur- 
face plasmons can also be detected using electron energy 
loss spectroscopy (EELS) [2S]. As analyzed in ref. [26J, 
EELS yields a signal related to the integral of the LDOS 
along the normal to the interface. Hence, the distance 
dependence of the LDOS is lost in EELS experiments. 
EELS has been successfully applied to study LDOS of 
graphene on SiC [12]. A remarkable result is the obser- 



vation of strong coupling between surface phonon polari- 
ton on SiC / air and surface plasmon in graphene and the 
demonstration of the tunability of the coupling by chang- 
ing an applied potential to control the chemical poten- 
tial. Yet, as pointed out, the distance dependence of the 
LDOS cannot be observed using EELS. 




FIG. 1: Sketch of a dielectric substrate covered with a 
graphene sheet and basic principle for the measurement of 
the LDOS above this system, obtained by deconvoluating the 
flux of electromagnetic field diffracted by an AFM tip toward 
a probe aperture. 

Here we focus in particular on the distance dependence 
of the LDOS. It had been shown [27] that the LDOS has 
dramatic spectral changes in the near field due to to the 
resonant contribution of the surface waves at the frequen- 
cies where e(w) + 1 = 0. These large spectral changes 
with distance in vacuum are a spectacular example of 
the non-invariance of spectra upon propagation in vac- 
uum first discussed in ref. [28]. In this paper, we show 
that when graphene is deposited on SiC the LDOS peak 
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positions strongly depend on distance. In addition, they 
also depend on chemical potential, paving the way to an 
active control of the LDOS. In particular we pay a specific 
attention to the description of the coupling mechanism 
between the graphene plasmon and the surface phonon- 
polariton supported by the SiC and highlight the finger- 
print of this coupling on the LDOS spectrum. 

The expression of the LDOS is usually decomposed in 
an electric and a magnetic contribution a,s p — pe + Pm, 
where 
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and pm is obtained from pE by exchanging the polar- 
izations TE and TM. In this expression k = ck/uj, k 
being the component of the wavevector parallel to the 
surface, p = \/\ — k? and p^{uS) = uP' /■k'^c^ is the LDOS 
in vacuum. The material properties of the surface are 
taken into account in the LDOS ([T|) through the Fresnel 
refiection coefficients tte and ttm for the two polariza- 
tions. In the general configuration of a substrate hav- 
ing a frequency-dependent dielectric permittivity e{uj) 
covered with graphene (having 2D surface conductivity 
c{uj)) these take the form ^1 [TU] 
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where p„i — \/ e(w) — equals c/w times the z com- 
ponent of the wavevector inside the medium. The ordi- 
nary case of a dielectric substrate is recovered by taking 
(j{uj) — 0, whereas e(a;) — 1 gives back the reflection 
coefficients of a suspended graphene sheet. 

Before calculating the LDOS, we need a model for the 
permittivity of SiC and for the conductivity 0(10) of 
graphene. For the permittivity we use a Lorentz model 
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6.7, UJL = 1-827 X 10i^rads"\ wt = 1.495 x 10"rads"^ 
and r = 0.9 X lO^^rads"^. This model predicts the ex- 
istence of a surface phonon-polariton resonance at fre- 
quency Wgpp ~ 1.787 X lO^^rads"^. As for the con- 
ductivity cr((u;), it can be written as a sum of an in- 
traband (Drude) and an interband contribution, respec- 
tively given by [lOj 



OB (w) 
cr/(w) : 



2 r 



■^^log(2cosh — 



e 
4/1 



«(t) 



-t-i- 



+00 



G(0 



(3) 



where G(x) = sinh(a;/fcBT)/[cosh(^/fc_Br) -|- 
cosh(a;/A:BT)]. The conductivity depends explicitly 
on the temperature T of the graphene sheet, for which 
we have chosen T = 300 K in our calculations. Moreover, 
it contains the chemical potential ^, which represents a 
varying parameter in our discussion, and the relaxation 
time r, for which we have used the value [TT] r = 10^^"^ s. 

We now turn to the analysis of the LDOS spectra. In 
Fig. [2]we show the electric, magnetic and total LDOS for 
a semi-infinite SiC surface (Fig. 2 a), for one layer of sus- 
pended graphene (Fig. 2 b) and for graphene on SiC (Fig. 2 
c). All the LDOS shown in Fig. [2] are calculated at a dis- 
tance of d = 50 nm from the surface, chosen to be smaller 
than the decay distance in air of SiC surface polaritons. 
The curves in Fig. |2|a) show the well-known peak [27] at 
w = Wspp due to the surface phonon polariton contribu- 
tion, almost entirely associated to the electric contribu- 
tion. The magnetic contribution, on the contrary, is at 
the origin of a secondary peak at w = cj^, approximately 
two orders of magnitude weaker, associated to the res- 
onance of the dielectric constant of SiC [30]. We now 
focus on a suspended sheet of graphene (Fig. |2jb)). In 
this case, we have observed that the electric part of the 
LDOS is manifestly dominating for any considered fre- 
quency: for this reason, we only represent the total LDOS 
both for suspended graphene and graphene-covered SiC. 
Note that the density of states is very broad as com- 
pared to the case of SiC and shows only one maximum 
less pronounced than for SiC. A remarkable property of 
graphene is the possibility of tuning the peak frequency 
by tuning the chemical potential /i. In particular, for 
the intermediate value \i = 0.2 eV, the peak matches the 
phonon-polariton SiC resonance. We are now interested 
in discussing what happens when graphene is deposited 
on SiC (Fig. |2jc)). For all considered /i, we observe the 
appearance of three peaks in the LDOS. One of them ap- 
pears at w = ujt- Its height and position do not depend 
on /i. This peak is due to s-polarized magnetic fields so 
that it is not associated with surface phonon polaritons 
and graphene plasmons. We now focus on the two other 
peaks. One of them keeps the memory of the SiC surface 
phonon-polariton resonance but is shifted as p increases 
and becomes broader. The other peak also shifts toward 
higher frequencies as p increases. 

We now discuss how in each considered configuration 
the LDOS results from the surface modes interacting 
within the graphene film. To this aim we study the poles 
of the reflection coefficients ([2]), by focusing in particu- 
lar on how the optical properties of graphene modify the 
shape of the surface modes of SiC. This phenomenon will 
be specifically discussed, as the LDOS presented in Fig. 
[2] as a function of p. It is well known [3T] that the opti- 
cal description of a polar material by means of a Lorentz 
model predicts the existence of surface modes only in TM 
polarization. For this reason we will limit our analysis to 
TM polarization also in the case of suspended graphene 
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FIG. 2: LDOS at a distance of d = 50 nm from a SiC substrate 
(a), a suspended graphene sheet (b) and a graphene-covered 
SiC substrate (c). In panel (a) the total LDOS (solid line) 
is decomposed in its electric (dot-dashed line) and magnetic 
contribution (dotted line). In panels (b) and (c), correspond- 
ing to suspended graphene and graphene-covered SiC respec- 
tively, only the total LDOS is represented, for three values of 
the graphene chemical potential: = OeV (blue dot-dashed 
line), \JL = 0.2 eV (red dashed line) and /i = 0.5 eV (black solid 
line) . 



and graphene on SiC. In order to derive the dispersion 
relation of surface modes in the three cases, we study the 
zeros of the denominator of ttm (see eq. (|2|). It is known 
that the dispersion relation depends on the choice made 
when searching the poles. As discussed in ref. |S2], when 
discussing LDOS, we look for poles with real wavevector 
k and complex frequency w = cj^ + *Wi • In this calculation 
scheme, provides the energy of the considered mode, 
while the inverse of uji defines its lifetime. We have per- 
formed this calculation for the three configurations in the 
same frequency region in which the LDOS has been dis- 
cussed, and for k varying between and 100 /im~^. This 
choice is associated to the fact that for a given distance 
d from the surface, the modes participating to the den- 



sity of states are smaller or of the order of [33] fcc — A , 
which for d = 50 nm gives k^ ~ 20 ^m~^ . 

The results are shown in Fig. [3] for the three config- 
urations and the same three choices of n used for the 
LDOS. The behavior of a Lorentz material in terms of 
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FIG. 3: Dispersion relations of surface modes for SiC (black 
dashed line), suspended graphene (green dotted line) and 
graphene-covered SiC (red solid line). The blue dot-dashed 
line represents the light cone in vacuum. The comparison 
is made for three different values of the graphene chemical 
potential (fi = 0(a), 0.2(b), 0.5(c) eV). 

surface-mode dispersion relation is well-known: it im- 
plies a branch following first (for small frequencies) the 
light cone, then approaching a horizontal asymptote at 
Lu — Wspp. For SiC, this asymptote is reasonably reached 
around k = 2//m~^. Thus, in the wavevector scale of 
Fig. [3] the SiC curve reduces as a matter of fact to a 
horizontal (dashed) line at w = i^spp- 

Let us now focus on suspended graphene. In this case, 
for any fi we observe the already known characteristic 
\/fc-like dispersion relation [TU]. Nevertheless, this ana- 
lytical dependence of w on A; is deduced by performing 
several approximations, and remarkably by taking into 
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account only the intraband contribution i7d(w) to the 
conductivity (see eq. ([3|). This simplification starts fail- 
ing when reducing the value of fi and increasing the value 
of w, even in the domains we are considering in this work. 
The result is a dispersion relation which has a similar 
shape, but increases always more slowly than the curve 
obtained for = CT£)(cl'). By repeating our calculation 
under this assumption and comparing the two results, we 
have observed that their relative difference reaches val- 
ues up to 115% for ^ = OeV and k = 100 /zm~^ (being 
instead 27% at fc = 20 /im"!), 7% for ^ = 0.5eV and 
k = lOO/im^i (1% at fc = 20^im-^). This shows that es- 
pecially for low values of the role played by interband 
transitions (described by ai{oj) in eq. (|3|) is essential 
for a precise quantitative calculation. Going back to the 
analysis of graphene curves, we see that they allow to 
readily explain the behavior of the LDOS shown in Fig. 
[2] b. First of all, the dispersion relation of plasmons on 
graphene does not have a flat asymptote so that we do 
not expect a peak in the LDOS. Yet, the LDOS displays 
broad peaks which move when varying fi. To explain 
their origin, we first observe that the integral describing 
the TM contribution to the electric density of states (eq. 
(jlj) is dominated by a factor exp(— 2kcliz/c) for large 
K as p is equivalent to iK. This clearly shows that there is 
a distance-dependent cutoff wavevector k^ = when 
observing at distance d. It is seen that the density of 
states has a peak value when exp{—2KU!d/c) is maxi- 
mum, i.e. for Kij/c — 1/d. This near-field filtering yields 
a d-dependent value of the most represented wavevector 
K. We can associate a frequency w to this peak k value 
using the dispersion relation. In summary, the frequency 
peak in the LDOS can be roughly estimated as the fre- 
quency associated to the largest participating wavevec- 
tor, as confirmed by comparing Figs. [2]and|3] 

We now finally discuss the case of graphene on SiC. 
The analysis of Fig. [3] shows the appearance of an anti- 
crossing, proving that we are in presence of a strong cou- 
pling between SiC phonon-polariton and graphene plas- 
mon as already shown |T5]. The interplay between the 
optical properties of the two materials gives rise to two 
separate branches of dispersion relation, whose proper- 
ties depend on the chemical potential. The most im- 
portant effect due to the SiC-graphene coupling is the 
disappearance of the horizontal asymptote at = (^spp- 
More specifically, for larger values of fi the high-frequency 
branch of the dispersion relation moves away from uj = 
Wspp for large values of k and has no longer a flat asymp- 
tote. Conversely, the lower branch increases and re- 
mains bounded by a horizontal asymptote uj = ujt- For 
/i = 0.2 eV and /i = 0.5 eV, the branch drifts toward 
higher values of frequencies, explaining the shift of the 
peak observed in Fig. [2] It is important to notice that 
the intersection between k — d^^ = 20 /im^^ and the 
dispersion curve provides a rough estimate of the posi- 
tion of the modified peaks. Once again, the position of 



the LDOS peak can be predicted from the dispersion re- 
lation, and the fact that this branch stays always below 
the one of suspended graphene corresponds to the fact 
that the peak in the LDOS is always at lower frequen- 
cies for graphene-covered SiC with respect to suspended 
graphene. 
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FIG. 4: LDOS for a SiC-graphene sample (/^ = 0.5 eV) for 
different values of the sample-tip distance d. The orange 
dot-dashed vertical line corresponds to a; = lot, whereas the 
dashed black vertical line is = ojspp. 

We have highlighted that the position of the frequency 
peaks in the LDOS can be predicted by the behavior 
of dispersion curves at the intersection points between 
k = d~^ and the different branches of the coupled modes. 
Hence, this analysis predicts d-dependent spectra of the 
LDOS. We display in Fig. |4]the LDOS for a SiC-graphene 
system with n = 0.5 eV for several distances d ranging 
from 1 nm to 10 fim. It is seen that the spectra varies dra- 
matically with distance. For SiC, the peak amplitude of 
the LDOS varies but its position remains fixed. Here, ow- 
ing to strong coupling between surface phonon polaritons 
and graphene plasmons, we observe peaks at frequencies 
which depend on the distance to the surface and on the 
chemical potential. This is a clear consequence of the in- 
terplay between strong coupling and near-field effects. To 
summarize, at a given distance d, we expect a dominant 
contribution to the LDOS of the mode with k ~ d^^. 
To this value correspond two frequency peaks given by 
Fig. |3] For instance, for the smallest considered distance 
d — Inm in Fig. [2j we have a cutoff kc ~ 1000 /im^^, 
the quasi-horizontal shape of the lower branch of coupled 
modes (see Fig. [sj^c)), gives rise to a pronounced peak in 
proximity of a; = ujt- In Fig- [4) we see that the LDOS 
spectrum displays two peaks due to strong coupling for 
d > 10 nm. On the contrary, beyond 1 fim the contribu- 
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tion of surface phonon polaritons and surface plasmons to 
the LDOS is reduced by five orders of magnitude so that 
the spectral peaks disappear. We conclude by emphasiz- 
ing that the distance dependence of the spectrum has to 
be accounted for when designing experimental measure- 
ments of the LDOS. 

In this paper, the combined analysis of the LDOS and 
dispersion relations has demonstrated a strong depen- 
dence of the spectrum of the LDOS on both the distance 
and the chemical potential. This shows that graphene- 
based heterosystems are tunable and actively control- 
lable materials [M]. We have also shown that the pres- 
ence of a graphene sheet significantly broadens the LDOS 
spectrum of the substrate. These LDOS properties offer 
new ways for controlling thermal management of nano- 
objects in lab-on-chip devices or modifying the sponta- 
neous emission rate of quantum emitters. 
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